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The calculational precision of e + e~ — > t + t~ and e + e~ — > fi + production cross sections in 
electron-positron annihilations at ^/s = 10.58 GeV is studied for the KKMC Monte Carlo simulation 
program, modified to include contributions from recent implementation of the hadronic part of 
vacuum polarization. We determine a(e + e~ — > t + t~) = (0.919±0.003) nb and a(e + e~ — » ^i + /i~) = 
(1.147 ± 0.005) nb, where the error represents the precision of the calculation. 
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I. INTRODUCTION 

At the present generation B-Factories lepton pair pro- 
duction, e + e~ — > t + t~ and e + e~ — > occurs at ap- 
proximately the same rate as the e + e~ — > BB process. 
Given the high statistics of the available data, precise 
knowledge of the r-pair and /x-pair production cross sec- 
tions in e + e~ annihilations at a center-of-mass energy y/s 
= 10.58 GeV is necessary for a number of precision mea- 
surements, particularly for measurements of the branch- 
ing fractions of r decays and luminosity determinations 
using the counting of e + e~ — ► events. The Belle 
experiment currently uses a r-pair cross section of 0.89 nb 
and quotes a 1.4% error on the number of r-pairs pro- 
duced for a given amount of integrated luminosity where 
the 1.4% is intended to account mainly for the luminos- 
ity measurement error [l], Q ■ The BABAR experiment also 
uses a 0.89 nb cross section but assigns a 2% error to this 
figure based on a simple comparison of the cross sections 
from the KKMC 3] and K0RALB Q Monte Carlo simula- 
tion programs [5j . Although the 2% figure is adequate in 
searches for rare or forbidden processes, such a large er- 
ror currently limits the precision of r branching fraction 
measurements at the B-factories 1 . 

In principle, the cross sections determined with KKMC 
are the most precise calculations currently available, but 
the cross sections and their precisions as determined with 
this code at energies far below the Z-pole have not been 
reported. Only tests for a/s > 40 GeV have been available 
in the literature 0, d, [t| and the code cautions users of 



An alternative approach that is insensitive to r-pair production 
cross section is to measure the r-branching ratios relative to the 
t electronic branching fraction which is measured to ~ 0.3% 
accuracy, as has been done in a recent publication Q. However, 
in such an approach, the corresponding Monte Carlo generators, 
e.g. TAU0LA Q, would need to be tested at the precision of 0.1% 
level in the presence of experimental cuts, and such a study is 
beyond the scope of the present paper. 



its potential lack of reliability for energies below y/s = 
40 GeV because of insufficient study in that regime. 

This unsatisfactory situation should be compared with 
predictions of KKMC for high energies, where a precision 
at the 0.1-0.4% level was achieved [H, Q. In this paper 
we report on straightforward studies of the [i and r pair 
production cross-section calculations of KKMC and K0RALB 
at ^fs = 10.58 GeV with a goal of providing, for the first 
time, a realistic estimate of the precision of the cross- 
section calculation at energies relevant to the B-factories 
relying on, but not repeating, the detailed high energy 
investigations. The cross sections are calculated using 
code that includes QED and electroweak radiative cor- 
rections to different orders in the fine structure constant 
a and using different approximations for vacuum polar- 
ization. K0RALB calculates the matrix element to order 
a including interference between initial (ISR) and final 
state radiation (FSR). KKMC calculates the cross section 
to order a 2 log^/m 2 ) and includes exponentiation and 
ISR-FSR interference. Various implementations of the 
vacuum polarization can be used in both K0RALB and 
KKMC programs. 

The results presented here rely on the extensive stud- 
ies [H ||| used to determine the precision of the calcu- 
lations in the KKMC as well as KORALZ [To| software. 
The conclusion of those studies is that cross sections 
calculated with KKMC at the Z-pole and LEP II ener- 
gies (y^s ~ 200 GeV) are at the precision level of 0.2%. 
Fortunately, the physics environment at the higher en- 
ergies where the detailed tests have been made is far 
more complex than at the T(45). For example, un- 
like at the T(4S I ), the higher energy tests had to ad- 
dress issues related to the radiative return to the Z- 
pole, including effects related to initial and final state 
radiation bremsstrahlung interference; scattering-angle 
dependencies, of particular interest because of the im- 
portance of forward-backward asymmetry measurements 
at the higher energies; and high-energy electroweak and 
QCD corrections. Therefore it is only necessary to verify 
that the technical precision determined in detail at the 
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higher energies for a wide class of effects can be reliably 
extended to the lower energies of the B-factories. This 
has been accomplished using relatively simple tests. In 
this way we achieve the goal of providing a realistic cross- 
section error without having to repeat the previous high 
energy studies. However, because of the different energy 
scales, the high energy studies of several effects, such as 
the vacuum polarization, are not reliably extended to the 
T(4S), and these have been examined separately. 

As it turns out, our current knowledge of the lepton 
pair production cross sections is in fact partially limited 
by the treatment of vacuum polarization in this calcula- 
tion. In this paper, we study their contributions using 
different parameterizations of the hadronic part of the 
vacuum polarization, as well as our knowledge of the ra- 
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Other sources 



of theoretical uncertainty arise from the implementation 
of initial and final state radiation including their inter- 
ference and electrowcak corrections, effects from virtual 
pair-productions of additional fcrmions, and the impact 
of low-energy resonances on the cross section. 

At the level of precision targeted in this work, cross 
sections for r-pair and /x-pair processes are adequately 
studied with 10 6 generated events. Because of the larger 
contribution from radiative returns to very low energy 
in /i-pair events, cross sections are also studied with an 
angular acceptance in the center-of-mass for both muons 
of [22°, 158°] and with effective center-of-mass energy, 
\fp , greater than 10% of 10.58 GeV which corresponds 
more closely to the actual experimental acceptance. 

Our paper is organized as follows. In Section [II] we dis- 
cuss radiative corrections to the total cross section due 
to hadronic vacuum polarization. We have found that 
changes to the public versions of the two programs are 
necessary in order to achieve the goal of obtaining a re- 
duced systematic error. Section HlTl is devoted to verifying 
that the high energy estimates of the systematic error of 
KKMC photonic (and acceptance dependent) corrections 
for the total cross section remain valid at our energies 
as well. Sections IIVI and [V] discuss the impact of in- 
terference and pair production corrections, respectively. 
Section IVT1 describes how we address the effects of vector 
resonances. Section fVIII discusses the error on the ratio 
of the t + t~ and /i + /i~ cross sections as calculated by 
KKMC. The estimations of the systematic error of theoret- 
ical predictions for the r-pair and /i-pair cross sections 
are summarized in Section IVIIII 



II. TREATMENT OF VACUUM 
POLARIZATION 

The particular vacuum polarization implementation is 
controlled in the KKMC code using the IHVP flag following 
the DIZET implementation of the calculation [ll[ , with 
improvements summarized in section 4.1.3 of Ref. [l2j 
and section 4.132 of Ref. [13j. The default modeling of 
vacuum polarization in KKMC corresponding to IHVP flag 



= 1 |l4[, suggested in Ref. [12j was not optimized for 
the low energy applications. In the present implemen- 
tation, the hadronic part of the vacuum polarization is 
studied using the experimental knowledge of the ratio 
R = a(e + e~ — ► qq)/a(e + e~ — ► /i + /i~) relevant for a 
precise determination of a(mz). The hadronic part of 
the vacuum polarization is not calculated in this option 
in DIZET for center-of-mass energies less than 40 GeV. 
Thus, at y/s= 10.58 GeV, this option gives a ~ 3% lower 
fcrmion-pair cross section than those obtained from other 
DIZET options that include the hadronic component 
of the vacuum polarization. Other available options in 
DIZET (IHVP flag = 2 [I| and = 3 [l6j]) use experimen- 



tal knowledge of the ratio R 



/i + fi~) measured more than 15 years ago. Improvements 
in the measurements of the ratio R from the Crystal Ball 
experiment (l7j in e + e~ annihilations at y/s between 5.0 
and 7.4 GeV, and more recent measurements by the BES 
experiment [18| at y/s between 2 and 5 GeV, are not 
included for these options. Improvements on the calcu- 
lation of the hadronic part of the vacuum polarization 
at y/s — niz obtained by including these new measure- 
ments of R have been documented recently [l9| . Follow- 
ing these developments, we have introduced a new option 
in DIZET implemented with IHVP flag = 4 to calculate 
the contribution of vacuum polarization using the routine 
"REPI" [2(| , where these improvements on the measure- 
ment of R have been incorporated. 

We evalute the uncertainties associated with using 
IHVP flag = 4. The routine REPI uses a simple 
parametrization of the hadronic contribution to the vac- 
uum polarization in the ^-channel and can be used safely 
in the s-channel outside the energy regime dominated 
by individual resonances. Special care must be taken 
to account for the large positive and negative fluctu- 
ations in the contribution to the hadronic part of the 
vacuum polarization in the s-channel across the reso- 
nance as illustrated in FIG. 1 of Ref. We have 
calculated the differences between the s-channel and t- 
channcl contribution of the hadronic part of the vacuum 
polarization to the r and /i-pair production cross sec- 
tion, weighted by the distribution of the invariant mass 
of the lepton-pair system (y/P) in bins of 10 MeV. Be- 
low and above the resonances, these differences are pos- 
itive and negative, respectively, and their effects tend to 
cancel when integrated over y/s'. At the T(45) these 
differences cancel partially, and the size of this cancela- 
tion is sensitive to the variation of the central value of 
the center-of-mass energy distribution of colliding beams 
for different periods of data taking as well as the collider 
beam energy spread. We assume a beam energy spread of 
4.6 MeV (RMS), which is characteristic of the center-of- 
mass ener gy s pread at BABAR experiment at the PEP II 
B-factory [22| ■ Based on the results of these calculations 
we assign an uncertainty of 0.18% on the r-pair cross 
section and 0.22% for the /i-pair cross section associated 
with the implementation of the vacuum polarization in 
KKMC with IHVP flag = 4 in DIZET. These uncertainties 
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are also valid for the Belle experiment, since the center- 
of-mass energy spread of the KEK B-factory is similar to 
that of the PEP II B-factory [H. 

In order to check the technical precision of this new 
implementation of vacuum polarization, we have also cal- 
culated the cross section with KKMC using different IHVP 
options in DIZET as well as with corresponding calcula- 
tions using KORALB with this new IHVP flag = 4 option 
implemented. The cross sections for e + e~ — > t + t~ and 
e + e~ — * calculated with KKMC are summarized in 

Table Q] along with their Monte Carlo statistical errors, 
which are at the level of 0.02% for r-pairs and 0.03% 
for the //-pairs. The cross sections for /i-pairs with the 
cuts mentioned in the previous section are denoted with 
"acc" and "cuts" subscripts, respectively, for the cases 
when only the angular acceptance cuts are applied and 
when both the acceptance and the Vs* > 0.1 y/s cuts are 
applied. For r-pair production, we report the total cross 
section. 

The T-pair cross sections from IHVP flags = 2, 3 and 
4 agree to within 0.16%. The /i-pair cross sections agree 
within 0.16% and 0.21% for the no-cuts case, cr(/i + /t~), 
and all cuts case, a cu t s (^ + respectively. As is ev- 
ident from the distribution of the lepton-pair invariant 
mass, a/s 7 , presented in FIG. [TJ the /x-pair cross section 
has a larger contribution from the vacuum polarization 
at low v / s'. Differences in the cross sections are expected 
because different measurements of R are used in the cal- 
culations of the vacuum polarization in Refs. (TH [l6l. l2C| . 
as described above. Therefore we do not assign any addi- 
tional uncertainty resulting from this cross check, when 
IHVP flag = 4 is used as the default option. 




Invariant Mass of fermion pair system (GeV) 



FIG. 1: Invariant Mass of r + r" 
IHVP flag = 4. 



and /t + /t pairs (y/s') with 



As an additional cross check on the introduction of 
the vacuum polarization into the code, we have stud- 
ied the technical robustness achieved in the implemen- 
tation of IHVP flag — 4 by comparing the cross sec- 



tions calculated with KORALB and KKMC programs using 
the different configurations given in Table HH Turning 
off the bremsstrahlung radiative corrections allows us to 
isolate the vacuum polarization corrections. In the case 
of KORALB, the vacuum polarization increases the Born- 
level cross section (ctborn LB ) by 27\Le(n 77 (s' = s)), where 
n 77 is the photon self-energy function for the vacuum po- 
larization correction. On the other hand, for KKMC, the 
vacuum polarization factor is |1 — n 77 (s')| 2 , where s' is 
only determined from ISR, so that for the case of FSR 
alone, s = s' . Therefore, we expect 
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This relation between KKMC and KORALB is found to be 
valid to better than the statistical precision of the Monte 
Carlo calculations for both r-pairs and //-pairs with and 
without additional cuts on acceptance and \fs' " . 

These results are consistent with our expectations. 
However, in general when ISR is turned on in KORALB 
s 1 = s is always used for the vacuum polarization calcu- 
lation, whereas, in KKMC s', more correctly, varies from 
event to event as expected in the presence of ISR. In 
the absence of ISR, the above relation comparing KKMC 
and KORALB cross sections is exactly valid, whereas in 
the presence of ISR, the difference in treatment of this 
effect between the two programs is small enough that the 
relation still provides a useful cross check. 



III. INITIAL AND FINAL STATE 
BREMSSTRAHLUNG 

Because bremsstrahlung is calculated to different or- 
ders in KKMC and KORALB, comparisons of the cross sec- 
tions calculated using the two software packages help 
study properties of the corrections, and as a conse- 
quence, establish the precision for the KKMC cross sec- 
tions. To that end, we compare results from Born-level 
calculations (KKMC or KORALB) , first order radiative effects 
(KORALB) and second order calculations with exponentia- 
tion (KKMC). For these studies, the electroweak and vac- 
uum polarization corrections are initially switched off. 

As a technical benchmark, we initially compare the 
Born-level e + e~ — > /i + /i~ cross sections calculated at 
xfs = 10.58 GeV using KKMC and KORALB when both elec- 
troweak and radiative corrections are switched off. For 
technical reasons [24j . the final state radiation in KKMC 
is not swiched off and therefore it is expected that the 

ratio ct^ rn /<£n LB is e q ual to 1 + (|)(f ) when no ex- 
perimental cuts are applied. To better than 0.01%, we 
find this to be the case. 

Having established that the Born-level cross sections 
are in good agreement, we proceed to use the KORALB 
cross sections with ISR-FSR switched on and with elec- 
troweak and vacuum polarization corrections switched off 
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TABLE I: r-pair and /i-pair cross sections with different IHVP flags using KKMC. Note: IHVP=1 does not inciude the hadronic 
part of the vacuum polarization. 



IHVP 


a(r + 0(pb) 


^(cV )(pb) 


o- acc (^ + /i~)(pb) 


a cu ts )(P D ) 


1 


891.45 ± 0.20 


1117.81 ± 0.35 


825.90 ± 0.26 


810.75 ± 0.25 


2 


919.85 ± 0.21 


1149.01 ± 0.37 


851.47 ± 0.27 


836.17 ± 0.27 


3 


920.17 ± 0.21 


1148.56 ± 0.37 


851.69 ± 0.27 


836.47 ± 0.27 


4 


918.66 ± 0.21 


1146.60 ± 0.36 


850.31 ± 0.27 


835.15 ± 0.27 



TABLE II: Code configurations. Note: for technical reasons, the KKMC code can only be run with final state radiation on. 





Software 


bremsstrahlung 


electroweak and vacuum 




Package 


configuration 


polarization corrections 


KORALB 

°"born 


KORALB 


ISR off, FSR off 


off 


KK 
°BORN 


KKMC 


ISR off, FSR on 


off 


^KORALB 


KORALB 


on with interference 


on 


KORALB 
^NO VP 


KORALB 


on with interference 


off 


KK 
°NO BREM 


KKMC 


ISR off, FSR on 


on 


KK 
^NO VP 


KKMC 


on with interference 


off 




KKMC 


on with interference 


on 


KK 
°"nO INT 


KKMC 


ISR/FSR on but interference off 


on 



« H v A P LB ). The quantity (l-<g BA P LB /<g RALB ) is a mea- 
sure of the contribution of first order bremsstrahlung to 
the cross section. The magnitude of the bremsstralung 
effects at the second and third orders in the exponential 
expansion can be estimated as (1 - Cvp IB /CT) 2 / 2 
and (1 - <g R v A p LB /cr^ LB ) 3 /6, respectively. For 
fT cuts (/i+/i-), o^g^WSSU™ = Ln 26 whereas for 
a(r+r-), <STp LB /<ohn LB = 1-1054. Taking 11% as 
the size of first-order correction from bremsstrahlung 
the second order in the exponential expansion series is 
~ 0.11 2 /2 = 0.0061. Therefore, since KKMC includes the 
second order, we expect a difference between the KKMC 
and KORALB cross sections to be ~ 1%. This 1% estimate 
is consistent with the value of 1.0072 for a™ VP /a^ R ^ B 
for cr cuis (/i + /i _ ), and somewhat larger than the value 
of 1.0012 for <j(t + t~), which also receives contributions 
from other effects, e.g. mass corrections. The next or- 
der of the leading log terms, which is not fully controlled 
in KKMC, can be estimated to be ~ 0.11 3 /6 = 0.0002. 
These tests confirm that KKMC and KORALB are yield- 
ing the expected relative cross section behaviour near 
y/s = 10.58 GeV. We can estimate the relative error as- 
sociated with the treatment of ISR and FSR to be equal 
to the size of the last term fully controlled in the approx- 
imation implemented in KKMC: a 2 \og(s/ml) = 0.0011. 
This next-to-leading-log term is larger than the third or- 
der leading log term, which is expected to be of order 
of ~ 0.11 3 /6, by about a factor of five. The above es- 
timates are rather naive as they assume an exponenti- 
ation pattern for all QED effects. Nevertheless, in our 
case the resulting relations hold rather well, pointing to 
a simple pattern of QED effects, and at the same time to 
the technical correctness of our calculations. This is to 
be compared with the 0.2% uncertainty assigned for the 
KKMC treatment of bremsstrahlung at LEP II energies @ , 



where further complications arise due to significant con- 
tributions from the radiative return to Z. 

Final state bremsstrahlung has mass-term dependen- 
cies that potentially are more significant for r-pair pro- 
duction. These would be expected to contribute a max- 
imum of — — — = 0.03% to the relative systematic error 
and therefore are negligible for this level of study and not 
considered further. 

As these tests assume KORALB is correctly calculating 
the matrix element to order a, it is useful to ensure that 
its behaviour associated with the technical parameter 
XK0, the minimal energy for the real bremsstrahlung pho- 
ton to be explicitly generated, is insensitive to the actual 
choice of this parameter [251 ]. To that end, we have ver- 
ified that reducing XPAR11=XK0 from its default value of 
0.01 by a factor of two changes total KORALB cross section 
by considerably less than 0.1%. 

Finally, as a technical cross check on the final re- 
sults, we also compare the KKMC cross-section calcula- 
tion with all radiative corrections turned on, <r KK , to 
a cross-section calculation assuming naive factorization 
of the bremsstrahlung and non-bremsstrahlung correc- 
tions: cr KK = <£ vp x (cr^ BREM /cr^ RN ). For ^-pairs 
with sfp > which is in the regime of very small 

radiation, we find that this factorization holds to better 
than the statistical precision of the Monte Carlo calcula- 
tions. 

These studies verify that the 0.1-0.2% error tag as- 
signed for photonic corrections in Refs. [1, Q for higher 
energies is also valid at y/s = 10.58 GeV. We adopt the 
higher, 0.20%, value as the uncertainty on the imple- 
mentation of ISR/FSR bremsstrahlung at the T(4S) and 
thereby avoid having to address the dependence on selec- 
tion cuts. 
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IV. INTERFERENCE EFFECTS 

Two kind of interference effects arc considered: the 
electroweak interference between 7 and Z* propagators, 
and the QED interference between intial and final state 
radiation arising from 77 box diagrams. The contribu- 
tion from the 7-Z* interference is smaller than the Monte 
Carlo statistical precision. Although the QED interfer- 
ence results in a forward-backward asymmetry of a few 
percent, the contribution to the total cross section is con- 
siderably smaller. The magnitude of the contribution of 
interference to the total cross section can be estimated 
from the ratio of the cross section with all corrections on, 
<7 KK , to the cross section calculated with all corrections 
on but with interference switched off, cr™ INT . For both 
^-pairs and r-pairs, 0" KK /cr££ INT = 1.0004 and is rela- 
tively insensitive to whether or not geometrically sym- 
metric cuts are applied. We conclude that the largest 
contribution to uncertainties from these interference ef- 
fects are of the order of 0.04%. 



V. PAIR-PRODUCTION AND VERTEX 
CORRECTION UNCERTAINTY 

The default calculations in KKMC do not include the 
effects of the emission of an extra fermion pair accom- 
panying the main process and the virtual fermion loop 
vertex correction. As these two effects are of opposite 
sign and comparable magnitude they will largely cancel 
each other and the error introduced by neglecting them 
both in the calculation is small. This large cancelation is 
a well-known feature of QED and has been numerically 
verified by studies at higher energies (l3| . 

To estimate the error this introduces on the total cross 
section, we activate the appropriate second order con- 
tribution to the vertex in KKMC as explained in Section 
4.54 of Ref. [l3j . and supplement the simulation with the 
additional sample of the four fermion final states, e.g. 
with the help of KORALW The effect of the vir- 
tual fermion loop vertex correction at the T(45), using 
appropriate options in KKMC, is 0.3%. However, because 
emission of real pairs largely cancels contribution from 
this vertex correction, we can conservatively assign half 
of this 0.3% as the error associated with ignoring both 
the pair-production and vertex virtual fermion loop cor- 
rections. We therefore assign a contribution of 0.15% to 
the uncertainty on the cross sections. 



VI. IMPACT OF RESONANCES 



production and decay of intermediate vector meson res- 
onances, J/#, #(25), Y(35), Y^S 1 ) etc, when accompa- 
nied by hard ISR: e.g. e + e~ — > J/#7 — > /i + /z _ 7. We es- 
timate the potential size of these contributions and add 
a systematic error that accounts for the fact that they 
are not included in the KKMC program. 

We calculate the cross section associated with the ra- 
diative return to each resonance below the T(45) using 
the narrow width approximation and a radiator function 
calculated to order a 2 [27L l28j . Updated resonance pa- 
rameters and leptonic branching ratios of reference [291 ] 
are employed. In the case of the J/#, we can cross check 
the calculation against the BABAR experiment's measure- 
ment of e+e~ — * J/#7 — > ^+^"7 [3(| and find they 
agree to 7%. Table Hill lists the vector meson resonances 
that give a small contribution to the cross section. Of 
the other vector meson resonances of potential interest 
(#(4415), #(4160), #(4039),#(3770), p(2150), p(1900), 
p(1700), 0(1680), w(1650), tti(1600), p(1450), w(1420), 
tti(1400), 0(1020), w(782) and p(770)) none contribute 
to the <7 cuts (/i + /i _ ) or <j(t + t~) cross sections at the lev- 
els of interest because they are either below the /z-pair 
y/P cut or r-pair threshold, or because they have a neg- 
ligible leptonic branching ratio. For the T(45), we use 
the measured cross section of 1.101 nb at the peak of the 
resonance [22f . The calculation of the contribution of the 
resonances on a C uts(n + £t _ ) employ the KKMC Monte Carlo 
simulation to estimate the impact of the cuts 2 . 

From this study it is evident that the vector meson 
resonances contribute at the 0.16% to the r-pair cross 
section and the 0.28% level to a cuts {^ + £*")• We assign 
these as conservative uncertainties on the KKMC cross sec- 
tion calculation because they are not included in the KKMC 
code. Note that if an s' cut is applied to a /i-pair selec- 
tion that removes the J/#(15) and #(25) resonances, 
e.g. v? > 0.4-^/s, then this component of the error on 
Cents (m + M~) is reduced to 0.15%, otherwise there will be 
an additional contribution of 0.12% from J/#(15) contri- 
bution to <J C uts{n + A* - )- If no cuts are applied, the error 
on a(n + fj,~) from this source is 0.33%. 



VII. UNCERTAINTY ON THE RATIO 

One of the means of determining the integrated lumi- 
nosity of a data set at the B-factories is to count /i-pair 
events. If this luminosity is used to establish the number 
of produced r-pair events in the sample, then the ratio 
(j(T + T~)/a cuts (p + fi~) appears in the denominator of the 
branching fraction calculation. In this case, the error on 



All of the studies so far discussed have assumed 
e + e~ — > r + r _ and e + e~ — > n + production to pro- 
ceed via 7 and Z* propagators. However, at y/s — 
10.58 GeV it is also possible for a non-negligible frac- 
tion of the final state lepton-pair yield to arise from the 



2 Note that the method employed to estimate the contributions 
to the resonances has been improved compared to that used in 
the first version of this paper (arXiv:0706.3235vl [hep-ph]) and 
results in a slightly larger determination of the total uncertainty. 
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TABLE III: Estimates of intermediate vector resonance contributions to the cross sections. 



Vector Resonance 


r total 

(MeV) 


BF( M +/0 
(%) 


Contribution to 

cr C uts{^ + >~) (%) 


BF(t+t-) 
(%) 


Contribution to 
a(r+r-) (%) 


T(4S) (10580) 


20.5 


0.0016 


< 0.01 


0.0016 


< 0.01 


T(3S)(10355) 


0.020 


2.18 


0.07 


2.29 


0.07 


T(2S) (10023) 


0.032 


1.93 


0.04 


2.00 


0.03 


T(1S)(9460) 


0.054 


2.48 


0.05 


2.60 


0.05 


*(2S)(3686) 


0.327 


0.74 


< 0.01 


0.30 


< 0.01 


J/#(1S)(3097) 


0.093 


5.93 


0.12 






Total 






0.28 




0.16 



this ratio contributes a theoretical systematic error asso- 
ciated with the overall normalization of the data sample. 
One can expect that some of the errors discussed in this 
paper will cancel in the ratio. 

For the vacuum polarization uncertainty, the only com- 
ponent of uncertainty that survives in the ratio is as- 
sociated with the contributions to ct(/j + /j~) below the 
t + t~ threshold because above the t + t~ threshold the 
vacuum polarization contributes to both r-pairs and fi- 
pairs in the same way, apart from small effects related 
to the exact description of the s' dependence above the 
threshold. We apply the same technique for evaluating 
the systematic error as described in Section ILT1 except we 
only consider the contributions up to the t + t~ thresh- 
old. This yields an uncertainty of less than 0.05%. The 
uncertainty of the determination of the hadronic vacuum 
polarization from the dispersion integral varies between 
0.01% and 0.02% below the r+r" threshold. These two 
uncertainties are independent and we assign the total un- 
certainty coming from the vacuum polarization in the ra- 
tio cr(r + T _ )/cr cnts (/i + ^ _ ) to be 0.05%. 

For the treatment of initial and final state 
bremsstrahlung, we conservatively take as the error 
on the ratio the same value of 0.2% assigned for the 
uncertainty on the absolute cross section. Almost 
certainly there are cancellations that further reduce this 
and which can be determined with additional study. 
Similarly, we assume the interference error on the ratio 
to be 0.04%. However, the pair-production and vertex 
correction uncertainties are 100% correlated and largely 
cancel in the ratio. Assuming the /i-pair selection 
employed for a luminosity determination applies a cut of 
\fs ! > 0.4-^/s, the error associated with the resonances 
will also cancel in the ratio, otherwise there will be an 



additional error of 0.12%. 

The total error we assign to the ratio 

a(T + T-)/a cuts (fi+n-) = 1.1000 is (0.05 © 0.20 © 
0.04 ©0.12)% = 0.24%. 



VIII. SUMMARY 

Adding the uncertainties from the effects discussed 
in this paper in quadrature we conclude that KKMC cal- 
culates the r-pair and /x-pair production cross sections 
at Vs = 10.58 GeV with a (0.18 © 0.20 © 0.04 © 
0.15 © 0.16)% = 0.35% and (0.22 © 0.20 © 0.04 © 0.15 © 
0.28)% = 0.44% relative calculational uncertainty, re- 
spectively. Using the IHVP flag = 4 option with the KKMC 
program as the central values, we obtain: cr(e + e _ — > 
t+t~) = (0.919±0.003)nb, a(e+e- -> = (1.147± 

0.005) nb, a cuts {e+e- = (0.835 ± 0.004) nb and 

a-(T + T-)/a cuts (fj, + iJ,-) = 1.100 ± 0.003. Given the ap- 
proach taken here to estimate the cross-section error, it 
is evident that, if required, the quoted precision can be 
further improved in the future by undertaking a more 
detailed, though well-defined, analysis. 
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